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Abstract. We study the dynamics of the FLRW flat cosmological models in which the vacuum energy density varies with 
time, A(t). In particular, we investigate the dynamical properties of a generalized vacuum model and we find that under certain 
circumstances the vacuum term in the radiation era varies as A(z) cc (1 + z) 4 , while in the matter era we have A(z) oc (1 + zf 
up to z — 3 and A(z) — A for z < 3. The confirmation of such a behavior would be of paramount importance because it could 
provide a solution to the cosmic coincidence problem as well as to the fine tuning problem, without changing the well known 
(from the concordance A-cosmology) Hubble expansion. 
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1. Introduction 

The analysis of the available high quality cosmological data 
(supernovae type la, CMB, galaxy clustering, etc.) have con- 
verged during the last decade towards a cosmic expansion his- 
tory that involves a spatial flat geometry and a recent accelerat- 
ing expansion of the universe (Spergel et al. 2007; Davis et al. 
2007; Kowalski et al. 2008; Komatsu et al. 2009 and references 
therein). This expansion has been attributed to an energy com- 
ponent (dark energy) with negative pressure which dominates 
the universe at late times and causes the observed accelerating 
expansion. The simplest type of dark energy corresponds to the 
cosmological constant (see for review Peebles & Ratra 2003). 
The so called concordance A model fits accurately the current 
observational data and thus it is an excellent candidate to be the 
model which describes the observed universe. 

However, the concordance model suffers from, among oth- 
ers (cf. Perivolaropoulos 2008), two fundamental problems: (a) 
the fine tuning problem ie., the fact that the observed value of 
the vacuum density (pa = Ac 2 /8ttG) is more than 120 orders 
of magnitude below that value found using quantum field the- 
ory (Weinberg 1989) and (b) the coincidence problem ie., the 
matter energy density and the vacuum energy density are of the 
same order prior to the present epoch, despite the fact that the 
former is a function of time but the latter not (Peebles & Ratra 
2003). Attempts to solve the coincidence problem have been 
presented in the literature (see Egan & Lineweaver 2008 and 
references therein), in which an easy way to overpass the coin- 
cidence problem is to replace the constant vacuum energy with 
a dark energy that evolves with time. The simplest approach 
is to consider a tracker scalar field <p in which it rolls down 
the potential energy V(<p) and therefore it could mimic the dark 
energy (see Ratra & Peebles 1988; Weinberg 1989; Turner & 



White 1997; Caldwell, Dave & Steinhardt 1998; Padmanabhan 
2003). Nevertheless, the latter consideration does not really 
solve the problem because the initial value of the dark energy 
still needs to be fine tuned (Padmanabhan 2003). Also, despite 
the fact that the current observations do not rule out the possi- 
bility of a dynamical dark energy (Tegmark et al. 2004), they 
strongly indicate that the dark energy equation of state param- 
eter w = PdeIpde is close to -1 (Spergel et al. 2007; Davis et 
al. 2007; Kowalski et al. 2008; Komatsu et al. 2009). 

Alternatively, more than two decades ago, Ozer & Taha 
(1987) proposed a different pattern in which a time varying A 
parameter could be a possible candidate to solve the two funda- 
mental cosmological puzzles (see also Bertolami 1986; Freese 
et al. 1987; Peebles & Ratra 1988; Carvalho, Lima & Waga 
1992; Overduin & Cooperstock 1998; Bertolami & Martins 
2000; Opher & Pellison 2004; Bauer 2005; Barrow & Clifton 
2006; Montenegro & Carneiro 2007 and references therein). 
In this cosmological paradigm, the dark energy equation of 
state parameter w is strictly equal to -1, but the vacuum en- 
ergy density (or A) is not a constant but varies with time. Of 
course, the weak point in this ideology is the unknown func- 
tional form of the A(f) parameter. Also, in the A(f) cosmologi- 
cal model there is a coupling between the time-dependent vac- 
uum and matter Wang & Meng 2005; Alcaniz & Lima 2005; 
Carneiro S. et al. 2008; Basilakos 2009; Basilakos, Plionis & 
Sola 2009). Indeed, using the combination of the conservation 
of the total energy with the variation of the vacuum energy, 
one can prove that the A(f) model provides either a particle 
production process or that the mass of the dark matter parti- 
cles increases (Basilakos 2009 and references therein). Despite 
the fact that most of the recent papers in dark energy studies 
are based on the assumption that the dark energy evolves in- 
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dependently of the dark matter, the unknown nature of both 
dark matter and dark energy implies that at the moment we 
can not exclude the possibility of interactions in the dark sector 
(eg., Zimdahl, Pavon & Chimento 2001; Amendola et al. 2003; 
Cai & Wang 2005; Binder & Kremer 2006; Das, Corasaniti, 
& Khoury 2006; Olivares, Atrio-Barandela & Pavon 2008 and 
references therein). 

The aim of this work is along the same lines, attempting to 
generalize the main cosmological properties of the traditional 
A-cosmology by introducing a time varying vacuum energy 
and specifically to investigate whether such models can yield 
a late accelerated phase of the cosmic expansion, without the 
need of the required, in the classical A-model, extreme fine tun- 
ing. The plan of the paper is as follows. The basic theoretical 
elements of the problem are presented in section 2, 3 and 4, by 
solving analytically [for a spatially fiat Friedmann-Lemaitre- 
Robertson-Walker (FLRW) geometry] the basic cosmological 
equations. Also in these sections we prove that the concordance 
A-cosmology is as a particular solution of the A(f) models. In 
section 5 we place constraints on the main parameters of our 
model by performing a likelihood analysis utilizing the recent 
Union08 Snla data (Kowalski et al. 2008). Also, in section 5 
we compare the different time varying vacuum models with the 
traditional A cosmology. In this section we also treat analyti- 
cally, the basic cosmological puzzles (the fine tuning and the 
cosmic coincidence problem) with the aid of the time varying 
A(f) parameter. Finally, in section 6 we draw our conclusions. 

2. The time dependent vacuum in the Expanding 
Universe 

In the context of a spatially flat FLRW geometry the basic cos- 
mological equations are: 



Ptot = Pf + Pa = 3H 
and 

d(p f +p A ) 



dt 



+ 3H(p f + P f + p A + P A )=0 , 



where p/ is the density of the "cosmic" fluid: 



Pf(t) = 



and 



p m (f) Matter era 
p r (t) Radiation era 



Pf(t)=Pp f 



Matter era fj = 



Radiation era /3 ■ 



1/3 



(1) 



(2) 



(3) 



(4) 



is the corresponding pressure. Also p A and Pa denote the 
density and the pressure of the vacuum component respec- 
tively. From a cosmological point of view, at an early enough 
epoch, the above generalized cosmic fluid behaves like radi- 
ation Pf = P r = p, /3 (J3 - 1/3), then behaves as matter 
Pf - Pm — (J3 — 0) and as long as Pa = -p\ it creates an 
accelerated phase of the cosmic expansion (see below). Notice, 
that in order to simplify our formalism we use geometrical units 
C&nG - c = 1) in which pa = A. In the present work, we would 
like to investigate the potential of a time varying A = A(f) 



parameter to account for the observed acceleration of the ex- 
pansion of the Universe. Within this framework it is interest- 
ing to mention that the equation of state takes the usual form 
of P A (f) = -pa (J) = -A(f) [see Ozer & Taha 1987; Peebles 
& Ratra 1988]. Also, introducing in the global dynamics the 
idea of the time-dependent vacuum, it is possible to explain the 
physical properties of the dark energy as well as the fine tuning 
and the coincidence problem respectively (see sections 5.1 and 
5.2). Using now eq.©, we have the following useful formula: 



Pf + 2>ifi + l)Hpf = -A 

and considering eq.([T]i we find: 



H + 



303+ 1) 



, A 
H 2 = - 
2 



(5) 



(6) 



where the over-dot denotes derivatives with respect to time. If 
the vacuum term is negligible, A(f) — > 0, then the solution of 
the above equation reduces to H(f) = 2(J3 + l)~ l /3t. Therefore, 
in the case of /3 = (matter era) we get the Einstein de-Sitter 
model as we should, H(t) = 2/3?, while for /3 — 1/3 we trace 
the radiation phase of the Universe ie., H(t) = 1/2/. On the 
other hand, if we consider the case of A(f) + then it becomes 
evident (see eq|5l) that there is a coupling between the time- 
dependent vacuum and matter (or radiation) component. 

Of course, in order to solve the above differential equation 
we need to define explicitly the functional form of the A(f) 
component. Note, that the traditional A = const cosmology 
can be described directly by the integration of the eq.© [for 
more details see section 3.1]. 

It is worth noting that the A(f) scenario has the caveat of 
its unknown exact functional form, which however is also the 
case for the vast majority of the dark energy models. Therefore, 
in the literature there have been different phenomenological 
parametrizations which treat the time-dependent Kit) function. 
In particular, Freese et al. (1987) considered that A(f) = 2>c\H 2 , 
with the constant c\ being the ratio of the vacuum to the 
sum of vacuum and matter density (see also Arcuri & Waga 
1994). Chen & Wu (1990) proposed a different ansatz in which 
A(f) oc a' 1 . 

Recently, many authors (see for example Ray, 
Mukhopadhyay & Meng 2007; Sil & Som 2008 and ref- 
erences therein) have investigated the global dynamical 
properties of the universe considering that the vacuum energy 
density decreases linearly either with the energy density or 
with the square Hubble parameter. Attempts to provide a 
theoretical explanation for the A(f) have also been presented 
in the literature (see Shapiro & Sola 2000; Babic et al. 2002; 
Grande et al. 2006; Sola 2008 and references therein). These 
authors found that a time dependent vacuum could arise from 
the renormalization group (RG) in quantum field theory. 
The corresponding solution for a running vacuum is found 
to be A(f) = co + c\H 2 (t) [where co and c\ are constants; 
Grande et al. 2006] and it can mimic the quintessence or 
phantom behavior and a smoothly transition between the two. 
Alternatively, Schutzahold (2002) used a different pattern 
in which the vacuum term is proportional to the Hubble 
parameter, A(a) oc H(a) [see also Carneiro et al. 2008], while 
Basilakos (2009) considered a power series form in H. It 
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is worth noting, that the linear pattern, A(a) oc H(a), has 
been motivated theoretically through a possible connection 
of cosmology with the QCD scale of strong interactions 
(Schutzhold 2002). In this context, it has also been proposed 
that the vacuum energy density can be defined from a possible 
link of dark energy with QCD and the topological structure of 
the universe (Urban & Zhitnitsky 2009). 

In this paper we have phenomenologically identified a func- 
tional form of A(a) for which we can solve the main differential 
equation (see eqj6j analytically. This is: 



A 



J2mt 



ym (t) = 3yH z (t) + 2mH(t) + 3n(J3 + 1 - y)e 

where the constants m and n are included for consistency of 
units (see below). Although, the above functional form was 
not motivated by some physical theory, but rather phenomeno- 
logically by the fact that it provides analytical solutions to the 
Friedmann equation, its exact form can be physically justified 
a posteriori within the framework of the previously mentioned 
theoretical models (see appendix A). 

Using now eq|7] the generalized Friedmann's equation (see 
eqO becomes 

fl = .30+l-y> g 2 +wg+ 3n(g+l 



2 2 - (8) 

and indeed, it is routine to perform the integration of eq.© to 
obtain (see appendix B): 

"303+ 1 -y) a/H 



H(t) = V^e ffl 'coth 



S(t) 



(9) 



(10) 



where 

f(eP--l)/» m*0 
\ t m = 

while the range of values for which the above integration is 
valid is n e (0, +oo) [for negative n values see the appendix]. 
Using now the definition of the Hubble parameter H = a I a, the 
scale factor of the universe a(t), evolves with time as 

"303+ 1 -y) V", 



a(t) = a\ sinh-^+'-T' 



S(t) 



(11) 



It is worth noting, that the relevant units of m + should cor- 
respond to time' 1 , which implies that m oc Hq. The parameter 
a i is the constant of integration given by 



\Pao/ 



(12) 



where pjo an d Pao are the corresponding densities at the present 
time [for which fl(fo) =1]. 

In this context, the density of the cosmic fluid evolves with 
time (see eqQ]i as: 

Pf (t) = 3H 2 (t)-A ym (t) (13) 
or 

p f (t) = 3(1 - y)H 2 (t) - 2mH(f) - 3n(J3 + 1 - y)e m ' . (14) 

In the following sections, we investigate thoroughly whether 
such a generalized vacuum component in an expanding 
Universe allows for a late accelerated phase of the Universe 
and under which circumstances such an approach provides a 
viable solution to the fine tuning problem as well as to the cos- 
mic coincidence problem. 



3. The matter+vacuum scenario 

In a matter+vacuum expanding universe (p/ = p m ), we attempt 
to investigate the correspondence of the A(f) pattern with the 
traditional A-cosmology in order to show the extent to which 
they compare. In particular, we will prove that the Hubble ex- 
pansion, provided by the current time-dependent vacuum, is a 
generalization to that of the traditional A cosmology. Note, that 
in the present formalism the matter era corresponds to f3 = 0. 



,q\ 3. 1 . The standard A-Cosmology 



Let us first investigate the solution for (y, m) = (0, 0). The vac- 
uum term eq.O of the problem becomes constant and is given 
by Aoo(a) = A = 3n. In this framework, the Hubble function 
(see eqj9jl is 



H A (t) = a/ - coth 



3 /A 
2V3 f 



Now, using the well know parametrization 

A = 3n = 3//^O a O a = 1 - Q m 

the scale factor of the universe is given by 

. ,j3H ^t\ 
a A (t) = a\ sinh ? 



1/3 



where (see eqlTZb 



<X\ — 



PmO 



1/3 



\Pao / \ Oa / 
The cosmic time is related with the scale factor as 



f A (a) = 



3VOa#o 




(15) 



(16) 



(17) 



(18) 



(19) 



Combining the above equations we can define the Hubble 
expansion as a function of the scale factor: 

H A (a) = H [Cl A + n m a- 3 ] l/2 . (20) 

In principle, Hq and O m are constrained by the recent WMAP 
data combined with the distance measurements from the type 
la supernovae (SNIa) and the Baryonic Acoustic Oscillations 
(BAOs) in the distribution of galaxies. Following the recent 
cosmological results of (Komatsu et al. 2009), we fix the cur- 
rent cosmological parameters as Hq = 70.5km/sec/Mpc and 
Q. m - I — Oa — 0.21. The current age of the universe (a = 1) is 
?oa - 13.77Gyr, while the inflection point takes place at 



tlA 



-sin 'h 



3 VO A #o 



fl/A 



o„ 



20a 



1/3 



(21) 



Therefore, we estimate t;\ 0.51foA and «/a — 0.56. 

Finally, due to the fact that the traditional A cosmology is 
a particular solution of the current time varying vacuum mod- 
els with (y, m) strictly equal to (0, 0), the constant value n is 
always defined by eq.dToTl. Therefore, throughout the paper all 
relevant cosmological quantities are parametrized according to 
n = A //ft. 
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3.2. "The general" A(f) Model 



3.3. "The modified" A Model 



In this section, we examine a more general class of vacuum 
models with (y, m) + (0, 0) (hereafter A ym model). The Hubble 
expansion and the corresponding evolution of the scale factor 
are (see eq|9]and eqfTTb 



Hit) = V^a H e""coth 
and 

a(t) = a\ sinh 3 "-'-' 



3(1 - y) yfO^Ho 
2m 



(e mf - 1) 



(22) 



3(l-y)VnA«o 



2m 



(e m ' - 1) 



or 



t(a) = —In 
m 



1 + 



2m . _ u I a 
— sin h — 

3(l-y)VQl//o \«i 



3(l-y)/2 



Obviously, if (y, m) — > (0, 0) [or e m ' - 1 » mf] then the A ym 
model tends to the traditional A cosmology, which implies that 
the latter be considered as particular solution of the general A ym 
model. Thus, this limit together with eq.(fT2l provide that 



(25) 



Taking the above expressions into account, the basic cosmolog- 
ical quantities as a function of the scale factor become 



H(a) = Ho [1 + 8(a)] [Q A + a„fl" 3(1 " y) ] 1/2 
and 

A ym (a) = 3yH 2 + 2mH + 3flgfl A (l - j)V + g(a)f 
where 



g(a) = 



2m 



sin : h 




'A fl 3(l-y)/2 



(26) 



(27) 



(28) 



It is worth noting that if we take (y, m) = (0, m) with m + 
(hereafter mild vacuum model or An m ), the corresponding 
Hubble flow becomes: 



H(a) = [l+g(a)]H A (d) 



(29) 



Therefore, as long as the function g(a) takes small values 
[g(a) <K 1], the Ao„, model has exactly the constant vacuum 
feature due to H(a) « H A {a). In this context, utilizing eq.dTTl) 
we simply have 



A 0m (a) = 2mH(a) + 3// 2 Q A [l + g(a)] 2 



(30) 



Finally, the fact that the vacuum term has units of time" 2 
implies that the vacuum term is proportional to H ( , or the con- 
stant m has to satisfy the following scaling relation: m oz Hq 
(see also section 2). Therefore, in the far future the condition 
m oc Ho + represents a super-accelerated expansion of the 
universe because a(t) oc exp( ^" H °£_), 



In this case we consider (y, m) = (y, 0) with y + (hereafter 
A y o model). From eq.© we can easily write the corresponding 
Hubble flow as a function of time 



(31) 



H(t) = V^a Hq coth 



Using now eqs.( [T0l ITTb . the scale factor of the universe a(f), 
evolves with time as 



(23) fl (0 = a \ sinh 311 



3(l-y)VHlH ( 



where 



i/3(i-r) 



(32) 



(33) 



Inverting eq.d32l we estimate the cosmic time: 



t(a) 



3(1 -y) y/H^H 



-sin 'h 



o 




A a 3(l-y)/2 



(34) 



The corresponding inflection point [a(tj) = 0] is found to be 
2 



3(l-y)V«Affo 



-sin 'h 



1 -3y 



or 



(l-3y)Q„ 



2Qa 



i/3(i~r) 



(35) 



(36) 



which implies that the condition for which an inflection point 
is present in the evolution of the scale factor is y < 1 /3. 

As expected, for y« 1 the above solution tends to the con- 
cordance model, a r o(f) — > «a(0- Now from eqs. (l3"n[32"l ). using 
the well known hyperbolic formula coth 2 jc -1 = 1/ sinh 2 x, we 
have after some algebra that: 



H(a) = H [Q. A + Q m a 



-3(i-r)]i/2 



(37) 



From this analysis, it becomes clear that the Hubble expansion 
predicted by the A y o model extents well that of the usual A cos- 
mology. To this end, utilizing eq.(l27b we can obtain the vacuum 
energy density 



A r0 (fl) = 3yH 2 (a) + 30 A // 2 (1 - y) 



(38) 



As we have previously mentioned in section 2, the above 
phenomenological functional form (see eql38l is motivated the- 
oretically by the renormalization group (RG) in the quantum 
field theory (Shapiro & Sola 2000; Babic et al. 2002; Sola 
2008). Moreover, recent studies (see Grande et al. 2006 and 
Grande, Pelinson & Sola 2009) find that this solution alleviates 
the cosmic coincidence problem (see section 5. 1). To conclude, 
it is worth noting that at late enough times (a » 1) the above 
solution asymptotically reaches the de-Sitter regime A ~ H 2 , 
as far as the global dynamics is concerned. 
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4. The radiation+vacuum scenario 

In this section, we consider a universe that is a spatially flat but 
contains both radiation and a time vacuum term. This crucial 
period in the cosmic history corresponds to B = 1/3. Therefore, 
for clarity reasons in the following sections we re-formulate our 
approach by using pj = p r and Pf = p r /3. These restrictions 
imply that 

PfO _ PrO_ _ Q_ 

Pao Pao 

where, Q,- =s 10~ 4 is the radiation density parameter at the 
present epoch derived by the CMB data (see Komatsu et al. 
2009). Within this context, based on eqs.©, dTTb and ( fT2l we 
present briefly the following cosmological situations: 



- radiation +constant vacuum: (y, m) 

factor is 

i 

a(t) = ^\\inh? (^H t) ■ 



(0,0): The scale 



(39) 



Owing to the fact that in this period t <K 1, the above so- 
lution reduces to the following simple analytic approxima- 
tion: 



a(t) * (2 VflrffoO 



a 1 

with H(t) ee - ~ — 
a It 



(40) 



- radiation+general vacuum: (y, m) + (0, 0): This general 
scenario provides 



a(t) = - sinh. 



2yi VQ^/Zo 



(e"" - 1) 



(41) 



where j\ = 1 - 3y/4. The vacuum component as a function 
of time (see eq|7]i is 



Aym(?) 



A ym (a) 



4(1 -yi) 
4y2f 2 



4( 1 - yi )Q r Hl 2m y/QHl 



(43) 



It is very interesting the fact that during the radiation epoch 
A ym (a) °c fiT 4yi . For small values of y or j\ ^ 0(1), the 
latter relation implies that as long as the scale factor tends 
to zero the vacuum term goes rapidly to infinity (see section 
6). In the case of (y, m) = (0, m) [or j\ - 1], the vacuum 
term (see eqsl42landl43l varies with time as 



(44) 



m 2m V^^o 
AomW * - ~ ~ 



Now the vacuum component evolves as A r o(a) °c gT 2 , in 
agreement with the Chen & Wu (1990) model. 
- radiation+modifled vacuum: (y, m) = (y, 0), y + 0: In 
this cosmological model we have 



(45) 




io- 6 10- 5 0.0001 0.001 0.01 0.1 1 
scale factor 



Fig. 1. Upper Panel: The evolution of the proximity parameter 
for the A r o cosmological model. Note, that the scale factor is 
normalized to unity at the present time. The lines correspond 
to y = 0.004 (solid) and y = 0.03 (dashed). Bottom Panel: The 
evolution of the radiation, matter and vacuum density consider- 
ing different kind of vacuums (after fitting the constants using 
the Union08 Snla data and Q m = 0.27, H = 70.5Km/s/Mpc). 
I) traditional A-cosmology: radiation density (open triangles), 
matter density (open squares) and constant vacuum density 
(open circles). II) modified A-cosmology, y + 0, A r o: radi- 
ation density (dashed line), matter density (long-dashed line) 
and vacuum density (solid line). Ill) The evolution of the mild 
vacuum, m + 0, Ao„, and IV) the evolution of the general vac- 



(42) uum, A ym (open stars). 



where yi = 1 — 3y/4. The approximate solution now be- 
comes 

a(t) * (2yj VfWof) I/2yi with H(t) * — . (46) 

2yif 

The vacuum component (see eqj7j evolves with time as 
4(1 -yi) 



Ay0(f) 



4yff 



2,2 



A(\-yi)a r Hl 
A r o(fl) ~ 4^ - A ym (a) 



(47) 



(48) 



Obviously, for a — > [yi - 0(1)] the vacuum energy 
density goes rapidly to infinity. 



5. Tackling the Cosmological puzzles 

As we have stated already in the introduction, there is a possi- 
bility for the vacuum energy to be a function of time rather than 
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Table 1. Numerical results. The l sf column indicates the vac- 
uum model used (the last two rows correspond to the fine tun- 
ing problem). Note, that the basic cosmological parameters 
were taken to be Q m = 0.27 and Ho = 70.5Km/sec/Mpc. 
Finally, the current age of the universe f has units of Gyr. 



Model 


7 


m/H a 


to 


A(tn) 


Ml pl ) 
A('n) 


A 








13.77 


1 


1 


A y0 


0.004 





13.82 


10 102 


10' 24 


Ao m 





2.4 x 10~ 3 


13.75 


10 51 


10 63 


Ay/jj 


0.004 


2.8 x 10~ 3 


13.80 


10 102 


10 124 



having a constant value. Therefore, in this section we compare 
the cosmic phases of the A(f) scenarios (described in the previ- 
ous sections) and the concordance A-cosmology. The aim here 
is to investigate the consequences of such a comparison on the 
basic cosmological puzzles namely cosmic coincidence prob- 
lem and fine tuning problem. 

5.1. The coincidence problem 

In order to investigate the coincidence problem we define the 
time-dependent proximity parameter of p m (a) (see eqfT4l> and 
Pa(ci) [see Egan & Lineweaver 2008 and references therein]: 



r(a) = min 



Pa(a) pm(a) 



Pm(a) p\(a) 



(49) 



Note, that in this work we use p\{a) = A(a) [see eqfT]. If the 
two densities differ by many orders of magnitude then r ^ 0. 
On the other hand if the two densities are equal the proxim- 
ity parameter is r — I. The current observational data shows 
that the proximity parameter at the present time (a = 1) is 
rn = = 7T- - 0.37. Therefore, a cosmological model 

may suffer from the so called coincidence problem if its prox- 
imity parameter is close to zero before the inflection point, 
r(a < a,) ~ 0. As an example, for the traditional A-cosmology 
we have r(a < 0.56) ~ 0. In contrast, if for a particular model 
we find that r(a < a/) = 0(1) then this model possibly does not 
suffer from the cosmic coincidence problem. 

In particular, suppose that we have a cosmological model 
which accommodates a late time accelerated expansion and 
it contains n-free parameters, described by the vector e = 
(e\,€2, e„). The main question that we should address here is 
'what is the range of input (ei, 62, e„) parameters for which 
the coincidence problem can be avoided?' Below we imple- 
ment the following tests. 

(i) We find the range of the free parameters of the consid- 
ered cosmological model that implies r ^ ro for at least two 
different epochs, one of which is precisely the present epoch. 

(ii) We know that for epochs between the inflection point 
and the present time a r < a < 1, the proximity parameter is 
r(a) > ro. As an example, for the traditional A-cosmology we 
have r(a) > 0.37. Thus, the goal here is to define the range 
of the free parameters in which at least a second region with 
r(a < ai) > ro occurs before the inflection point (a < aj). 

(ii) Once, steps (i) and (ii) are accomplished, we finally 
check whether the remaining parameters fit the recent Snla 



data, by performing a standard x 1 minimization. In this work, 
we use the so called Union08 sample of 307 supernovae of 
Kowalski et al. (2008). In particular, the x 1 function can be 
written as: 



307 



X 



./=! 



fi th (aj,e)-^(aj) 



(50) 



where a } ■, = ( 1 +z/)~ i is the observed scale factor of the universe, 
Zj is the observed redshift, [i is the distance modulus n = m- 
M = 51ogo?L + 25 and d\Sa, e) is the luminosity distance, given 

by 



Hoa X 



dx 



x 2 E(x) 



(51) 



where e is the vector containing the unknown free parameters 
and c is the speed of light (= 1 here). 

A cosmological model for which the present tests are suc- 
cessfully passed should not suffer of the coincidence problem. 
Below, we apply our tests to the current A(f) cosmological 
models (see also Table 1). 

- The modified vacuum model with e = (y, 0, ...0): We sam- 
ple the unknown y parameter as follows: y e (-1,1/3) 
in steps of 10~ 4 . We confirm that in the range of y e 
[0.004,0.03] the A y0 modefl satisfies both the criterion (i) 
and (ii) respectively. Also, we verify that this range of val- 
ues fits very well the Snla data, ^,„/dof 1.01. Notice, 
that for y > 0.03 the criterion (i) is not satisfied. As an ex- 
ample, in the upper panel of figure 1 we present the evo- 
lution of the proximity parameter for y = 0.004 (solid 
line) and 0.03 (dashed line). It is worth noting, that for 
0.1 < a < 0.34 (or 2 < z < 10) the vacuum density is 
low enough (r ~ 0) to allow galaxies and galaxy clusters 
to form (Garriga, Livio & Vilenkin 1999; Basilakos et al. 
2009). From now on, we will utilize y ^ 0.004 that cor- 
responds to the best fit parameter. To this end it becomes 
clear that the A y o model passes the above criteria and thus 
it does not suffer form the cosmic coincidence problem. 

- The mild vacuum model with e = (0, m, ...0): In this cos- 
mological model, we find that for m > 0.17//q, the corre- 
sponding age of the universe is to < 12.7Gyr. The latter 
appears to be ruled out by the ages of the oldest known 
globular clusters (Krauss 2003; Hansen et al. 2004). Using 
this constrain the unknown m parameter has an upper limit 
of 0.17//o and thus we perform the following sampling: 
m € [5 x 10 - 4 Ho,0.l7H ) in steps of 5 x IO^Hq. Within 
this range, we find that the required (i) and (ii) criteria are 
not satisfied. Thus, the Ao,„ cosmological model suffers of 
the coincidence problem. The resulting minimization pro- 
vides: m = 2A + _\ x m^Ho with ;tf m „/ dof - Note that 



1 Note, that from a theoretical viewpoint the predicted value of the 
y parameter is |y| = -p^ ^r, where M P is the Planck mass and M is an 
effective mass parameter representing the average mass of the heavy 
particles of the Grand Unified Theory (GUT) near the Planck scale, 
after taking into account their multiplicities. In the case of M ~ M P 
we can derive an upper limit of |y| < 1/127T (for more details see 
Basilakos et al. 2009). 
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the errors of the fitted parameters represent lcr uncertain- 
ties. 

The general vacuum model with e = (y, m, ...0): This 
vacuum cosmological model contains 2 free parameters. 
Using the previous mentioned sampling, we obtain that 
our main criteria for the A 



yirl 



scenario are full-filled for 

>-3i 



y e [0.004,0.02], m e [1.4 x 10~ 3 H ,9 x 10~ 3 H ] with 
X min /dof 6 [1.01,1.02]. Throughout the rest of the pa- 
per we will use the best fit parameters. These are: m - 
2.8 x 10- 3 H Q andy a 0.004 

In addition to the Snla data, we further check our statistical 
results using the dimensionless distance to the surface of the 
last scattering R = 1.71 ± 0.019 (Komatsu et al. 2009), and 
the baryon acoustic oscillation (BAO) distance at z = 0.35, 
A = 0.469 ± 0.017 (Eisenstein et al. 2005; Padmanabhan, et al. 
2007). We find that the above results remain unaltered. 



5.2. The cosmic evolution - fine tuning problem 

Using now our best fit parameters for the different kind of vac- 
uums, we present in figure 1 the corresponding normalized en- 
ergy densities, vacuum A(a)///?, matter p m (a)/H^ and radia- 
tion p r (a)/Hy as a function of the scale factor. We verify that 
both the A y o (solid line) and A ym (open stars) solutions are 
models that provide large values for the vacuum energy den- 
sity at early epochs, in contrast with the usual A cosmology 
(open circles), in which the vacuum energy density remains 
constant everywhere. Also, within a Hubble time (0 < a < 1) 
and for each (y, m) pair, we find the well known cosmic behav- 
ior for the matter density, p m (a) oc gt 3 and the radiation density, 
Pr(a) oc cT A respectively. As an example, in figure 1 we present 
the density evolution of the cosmic fluid for the A r o cosmologi- 
cal model: matter (dashed line) and radiation (dot-dashed line). 
For comparison we also plot the predictions of the traditional 
A cosmology: matter (open squares) and radiation (open trian- 
gles). From figure 1, it becomes clear that the radiation-matter 
equality takes place close to a rm a 3.7 x 10~ 4 a Sl r /£l m . For 
those vacuum models where m + (Ao m and A ym ), we verify 
that the behavior of their cosmic fluid (matter +radiation) devi- 
ates from the A solution in the far future (f » f ), since the 
exponential term e m ' in eq.(fT4li plays an important role in the 
global dynamics (see section 3.4 and below). 

In particular, for the A y o vacuum scenario (the same behav- 
ior holds for A ym ) we have revealed the following phases: (a) 
at early enough times (a < a rm ) the scale factor of the uni- 
verse tends to its minimum value, a — > 0, which means that 
the vacuum energy density initially goes quickly to infinity. So, 
as long as the scale factor increases the vacuum energy rolls 
down rapidly as A y o(a) oc a^ n [where j\ ~ O(l)]. This evolu- 
tion may treat the fine tuning problem. Indeed, for y e (0, 1 /3), 
we find that prior to the inflation point (f,„y ~ 10~ 32 sec), 
the vacuum energy density divided by its present value is 
A(f/„/)/A(?o) ~ 10 102 Finally, if we consider that the func- 
tional form of A(a) oc a~ 4yi is still valid during the Planck time 
{t p i ~ 10~ 43 sec), then A(f p/ )/ A(f ) ~ 10 124 (see the last rows in 
Table 1), and (b) in the matter era the vacuum density continues 
to roll down but with a different power law A y o(a) oc cr 3(1 ~ y) 
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Fig. 2. Upper Panel: Comparison of the scale factor provided 
by our A y o model with the traditional A cosmology (open 
points). Note, that we use Q. m = 0.27 and Hq = 70.5Km/s/Mpc 
model. In the bottom panel we present the deviation of the scale 
factors between the A y0 and A y ,„ model respectively. Note, that 
the scale factor is normalized to unity at the present time. 



and it tends to a constant value close to a ~ 0.25 (z ~ 3). 
Finally, for a > 0.25 the vacuum energy density is effectively 
frozen to its nominal value, A y o(a) a A = 3Qa//q, which im- 
plies that the considered time varying vacuum model explains 
why the matter energy density and the dark energy density are 
of the same order prior to the present epoch. The moment of 
radiation-vacuum equality occurs at a rv a 0.1 a {Q. r IQ.\) 1 ^ ■ 
Similarly, the moment of matter-vacuum equality takes place at 
a mv a 0.72 a (Q m /Q A )'/ 3 . From the observational viewpoint, 
in order to investigate whether the vacuum energy density fol- 
lows the above evolution, we need a robust cosmological probe 
at redshifts z > 3. In a recent paper (Basilakos et al. 2009), 
we have investigated how realistic it would be to detect differ- 
ences among the vacuum models. In particular, we have found 
that the Sunayev-Zeldovich cluster number-counts (as expected 
from the survey of the South Pole Telescope, Staniszewski 
et al. 2009, and the Atacama Cosmology Telescope, Hincks 
et al. 2009) indicate that we maybe able to detect significant 
differences among the vacuum models in the redshift range 
2.5 < z < 3 at a level of ~ 6 - 12%, which translates in number 
count differences, over the whole sky, of ~ 100 clusters (see 
figure 6 in Basilakos et al. 2009). 

Finally, in figure 1 we also show the evolution of the mild 
vacuum model Ao m (a) (dot line), in which y = 0. Briefly, 
we get the following dependence: (a) Ao,„ oc a~ 2y ' for a < 
a rm , while we estimate that Aom(^n/)/Ao m (^o) ~ 10 51 and 
Aom(^z)/Aom(*o) ~ 10 63 , (b) between a rm < a < 0.08 we have 
Aom x a~ 3 ^ 2 and (c) for a > 0.08 the Ao m becomes constant. 

We would like to end this section with a discussion on the 
evolution of the scale factor. In particular, our approach pro- 
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vides an evolution of the scale factor in the A y o model seen in 
the upper panel of figure 2 as the solid line, which mimics the 
corresponding scale factor of the A cosmological model (open 
points), despite the fact that they describe differently the vac- 
uum term. On the other hand, in the bottom panel of figure 2 
we present the corresponding deviation [(a ym - Oyo)/O)o]%, of 
the growth factors. It becomes evident, that within the range 
< Hot < 5 the evolution of the scale factor provided by the 
A ym model closely resembles, the corresponding scale factor of 
the A y o model (the same result holds also for the A cosmology). 
However, for models where m + 0, the situation is somewhat 
different in the far future. Indeed, for Hot > 5 the A ym (or An m ) 
cosmological scenario deviates from the A y o (or A) model by 
~ 5 - 10%. Thus, we conclude that the models with m + give 
a super-accelerated expansion of the universe in the far future 
with respect to those vacuum models where ra = 0. 

6. Conclusions 

The reason for which a cosmological constant leads to a late 
cosmic acceleration is because it introduces in Friedmann's 
equation a component which has an equation of state with neg- 
ative pressure, Pa = -pA- In the last decade the so called con- 
cordance A-cosmology is considered to be the model which 
describes the cosmological properties of the observed uni- 
verse because it fits accurately the current observational data. 
However, the traditional A cosmology suffers from two funda- 
mental puzzles. These are the fine tuning and the cosmic coin- 
cidence problems. An avenue through which the above cosmo- 
logical problems could be solved is via the time varying vac- 
uum energy which has the same equation of state as the tradi- 
tional A-cosmology. 

We wish to spell out clearly which are the basic assump- 
tions and conclusions of our analysis. 

- We are assuming a time varying vacuum pattern in which 
the specific functional form is: A(f) = 3yH 2 (t) + 2mH(t) + 
3n(J3 + 1 - y)e 2m/ , where (3 = (matter era) or = 1/3 
(radiation era), n = 3Qa//q, while the pair (y, m) charac- 
terizes the different types of vacuum. It is worth noting, 
that the above functional form includes the effect of the 
quantum field theory (for m = 0) [Shapiro & Sola 2000; 
Babic et al. 2002; Grande et al. 2006; Sola 2008] and it 
also extents recent studies (see for example Ray et al. 2007; 
Carneiro et al. 2008; Sil & Som 2008; Basilakos 2009). 
Notice, that we can easily prove that the cosmological con- 
stant is a particular solution of the general vacuum, that 
(y, m) = (0,0). Also we have investigated the following 
models: (a) modified vacuum in which (y, m) = (y, 0), mild 
vacuum with (y, m) = (0, m) and general vacuum in which 
(y, m) + (0, 0). In this framework, we find that the time evo- 
lution of the basic cosmological functions (scale factor and 
Hubble flow) are described in terms of hyperbolic functions 
which can accommodate a late time accelerated expansion, 
equivalent to the standard A model. 

- We find that that within the framework of either the mod- 
ified or general vacuum models the corresponding vac- 
uum term in the radiation era varies as A(a) oc aT A while 



in the matter dominated era we have A(a) ~ a 3 up to 
z = a - x - 1 3 while A(a) * A = 3Q A // 2 for z < 3. 
This vacuum mechanism simultaneously sets (a) the value 
of A at the present time to its observed value and (b) at 
the Planck time to a value which is 10 124 its present value 
[A(f p ;)/A(f ) ~ 10 124 ]. Additionally, we verify that our 
models appear to overcome the cosmic coincidence prob- 
lem. Finally, in order to confirm the above results, we need 
to define a robust cosmological probe at high redshifts 
(z > 3). In Basilakos et al. (2009) we propose that the future 
cluster surveys based on the Sunayev-Zeldovich detection 
method give some hope to distinguish the closely resem- 
bling vacuum models at high redshifts. 

Appendix 

In this appendix we provide a physical justification of the func- 
tional form of A(a) used in our paper. As we have already men- 
tioned in the section 2, the vacuum energy density can take 
several forms, depending on the theoretical approach. Briefly, 
the renormalization group from the quantum field theory intro- 
duces only even powers of H out of which the H 2 is the leading 
term (Grande et al. 2006; Sola 2008 and references therein). In 
another vein, the aforementioned possibility that the vacuum 
energy could be evolving linearly with H has been motivated 
theoretically through a possible connection of cosmology with 
the QCD scale of strong interactions (see Schutzhold 2002; 
Carneiro et al. 2008). In this framework, it has also been pro- 
posed a possible link of dark energy with QCD and the topolog- 
ical structure of the universe (Urban & Zhitnitsky 2009). The 
simplest approach therefore to introduce the effects of the DE is 
to consider a potential V(<p) Vo + m 2 (f> 2 /2, where the homoge- 
neous scalar field (f> obeys the Klein-Gordon equation. It is well 
known that for H const the corresponding (f> evolves with 
time as (p(t) (poe mt (where in general m is a complex number). 
In this context, one would expect that the functional form of the 
A(f) should contains also an additional term of (p 2 (t) oc e 2mf in 
order to take into account the possible link between dark energy 
and QCD. 

All the above options have merits and demerits. In the cur- 
rent paper, the functional form of A(f) is motivated by a com- 
bination of the above possibilities namely H 2 (t) [RG], H(f) 
[QCD] and e 2 '"' (dark energy). In particular, the linear com- 
bination reads as follows: 

A(f) = n x H 2 (i) + n 2 H(t) + n 3 e 2 '"' 

which obviously is very similar to the original (phenomeno- 
logically selected) form of A(f) (equation 7). Finally, from a 
mathematical point of view we can select the constants n\, ni 
and «3 to match with those presented in the original equation 7. 

Appendix 

With the aid of the differential equation theory we present so- 
lutions that are relevant to our eq.©. If one is able to have a 
Riccati differential equation which is given by the following 
special form 

^ = f(x)y 2 (x) + my(x) - ne 2m *f(x) (52) 
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then the general solution of eq.(|52l> for n > is 
y(x) = Vne""coth - \Jn f e mu f(u)du . 

J. c 

On the other hand, if n < then the solution of eq.(|52b is 



y(x) 



ne"~cot 



■Vw r 



e m "f(u)du 



(53) 



(54) 



Note, that in our formulation the function /(jc) is a constant: 
/(jt) = -306 + 1 - y)/2. Also, n < implies that Q„, > 1 (or 
A < 0). 
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